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Introduction 




The correspondence between M = 4 super symmetric Yanj 


5-Mills theory and string theory 


on 


an AdS$ x S5 geometry is the most studied example 


of a gauge/string duality [1-4]. 



After the impressive progress that has been made in computing anomalous dimensions of 
gauge invariant operators, see e.g. [5-13] and references therein, the focus has recently 
shifted to the analysis of planar scattering amplitudes. 

On the gauge theory side, the discussion was stimulated by the intriguing BDS formula 
of Bern, Dixon and Smirnov [14]. It encapsulates the known infrared and collinear behavior 
of n-particle maximally helicity violating (MHV) amplitudes in the planar approximation. 
The authors of [14] conjectured the BDS formula to determine the amplitudes at each loop 
order L > 2 up to some additive finite function of the kinematic variables. has 
been dubbed the 'remainder function'. Both string and gauge theory arguments suggest 
that the remainder function R^ vanishes for n = 4, 5. In the weakly coupled theory, 
perturbative computations have uncovered an important new symmetry of scattering am- 
plitudes [15]. Assuming this so-called dual conformal symmetry, the remainder functions 
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R(n) 

can only depend on conformal cross ratios, i.e. on conformal invariant combinations of 
the usual kinematic variables. Since there are no such cross ratios for n = 4,5, the corre- 
sponding remainder functions have to be trivial. In other words, dual conformal invariance 
predicts that the BDS formula is exact for n = 4, 5 to all loop orders. This prediction was 
confirmed by a string theory computation of the leading term at strong coupling. We shall 
say a bit more about the string theoretic analysis below. 

On the other hand, the remainder function R^ is known to be non-zero for n > 5 
and beyond one loop [16, 17]. A first explicit disagreement between the BDS formula 
and calculations of Yang-Mills scattering amplitudes in the leading logarithmic high energy 
approximation was pointed out in [18, 19]. This outcome has been confirmed by a numerical 
analysis of existing two- loop calculations for the six-point function [20]. 

The analysis in [18, 19] illustrates the importance of the high energy limit (Regge limit) 
which, in particular, exhibits much of the subtle analytic structure of multiparticle ampli- 
tudes. Let us briefly recapitulate the main results. There are three different kinematical 
regions which are of interest. In the usual 'physical region' all energies of the process are 
positive. The opposite case in which all energies are negative is referred to as the 'Eu- 
clidean region'. Most of the analysis we shall discuss in this work, however, deals with a 
third region that contains both positive and negative energies, see eq. (2.8) and Fig. 1 for a 
precise definition of this 'mixed region'. It is the mixed region that offers the clearest view 
on the discrepancy between the BDS formula and the existing high energy calculations of 
gluon amplitudes. In fact, the six-gluon scattering amplitude in the mixed region includes a 
Regge-cut contribution which cannot be reproduced from the BDS formula. For the physi- 
cal regime, the failure of the BDS formula manifests itself in a violation of the Steinmann 
relation. No disagreement is visible in the Euclidean region. 

These findings imply that the remainder function vanishes in the high energy limit 
of both the Euclidean and the physical region. Upon analytic continuation into the mixed 
region, however, the high energy limit should be nonzero and, in the leading logarithmic 
approximation, must coincide with the contribution given in [19]. This emphasizes the 
importance of the correct analytic structure. In a recent study [21], a convenient analytic 
expression [22] for the exact two- loop calculation of the six-point function [23, 24] was used 
to perform the relevant analytic continuation. The results are in full agreement with [19]. 
This settles the remainder function in the two-loop approximation, and it supports the 
all-order leading log generalization in [18, 19]. 

The leading logarithmic approximation in the mixed region exhibits an important fea- 
ture. As a function of the energy variable S2 (see Fig.l) the six-gluon amplitude Amp' 1 
contains a Regge cut term with a power-like behavior 

Amp' ~ , (1.1) 

where the (negative-valued) exponent E2 is the lowest eigenvalue of the BFKL color-octet 
Hamiltonian H2 for a two-gluon system. This term is missing in the BDS-amplitude and 

1 The / is placed on this amplitude to remind us the high energy limit is taken in the mixed rather than 
the physical region. 



-2- 



therefore should be connected with the remainder function . The BFKL color-octet 
Hamiltonian coincides with the Hamiltonian of an integrable open spin chain [25]. Similar 
results hold for scattering amplitudes involving n > 6 gluons. In a processes 2 — > 2n there 
exist Regge cut contributions with up to n gluons and the corresponding spin chains consist 
of n spins. Hence, the weak coupled theory provides direct evidence for integrability in the 
high energy behavior of planar scattering amplitudes. 

Having reviewed all these results from gauge theory it is natural to ask what string 
theory has to say about the high energy limit of the remainder function In order to 

understand how the issue can be addressed, we need to briefly sketch the development that 
was initiated by the work [26] of Alday and Maldacena. The main insight of this paper was to 
identify the leading strong coupling contribution to an n-gluon amplitude with the area A n 
of some 2-dimensional surface S n inside AdS§ . According to the prescription of [26] , S n ends 
on a piecewise light-like polygon in the boundary of AdS§. The light-like segments of this 
polygon are given by momenta pj of the external gluons. For n = 4 it is possible to find the 
surface explicitly and the resulting amplitude matches the prediction from the BDS formula. 
Constructing S n for n > 5, however, turned out to be a rather difficult problem, at least 
for finite n and generic choice of the external momenta. The issue was resolved through a 
series of papers [27-29] in which the area of S n is related to the free energy of some auxiliary 
quantum integrable system. More precisely, it was argued that A n may be computed from a 
family of functions Y s j with s = 1,2,3 and j = 1, . . . , n— 5. The latter can be determined by 
solving a set of coupled non-linear integral equations. Very similar mathematical structures 
are actually familiar from the study of ground states in 1-dimensional quantum integrable 
systems. Moreover, the functional A n resembles expressions for the free energy of such 
systems. 

We employ this auxiliary quantum integrable system to study the Regge limit of the 
remainder function at strong coupling. The non-linear integral equations for the Y- 
functions we have mentioned in the previous paragraph provide an explicit set of equations 
which enable us to take the high energy limit and to investigate the analytic continuation 
between the physical and the mixed region. As the main result, we determine the value of 
the exponent E2 at strong coupling, 



In passing to the mixed regime, we need to analytically continue the parameters of the 
auxiliary quantum integrable system. This continuation leads to a new system of equations 
with additional terms. The structure of the new integral equations and the corresponding 
modifications of the relevant energy functional are again familiar from the theory of quantum 
integrable systems. In that context they are used to describe excitations of the system. Our 
formula for the value of E2 at strong coupling assumes that we can perform the analytic 
continuation into the mixed regime within the parameter-space of the quantum integrable 
system and that the corresponding amplitude is given by the exponential of the energy 
functional, as in the Euclidean regime. 





-3- 



This paper is organized as follows. In section 2 we define the kinematic regions and 
the path of analytic continuation for the kinematic variables. Section 3 starts with a brief 
summary of relevant results from the work of Alday et. al. As a first exercise we then 
verify that the Regge limit of the remainder function vanishes in the physical regime 
(up to an irrelevant constant). Sections 4 and 5 contain the main part of our analysis. In 
these two sections we perform the analytic continuation, discuss the origin of the additional 
terms and calculate their contribution to the amplitude. For pedagogical reasons, we first 
illustrate the general arguments with a simpler example that is taken from the work [30, 31] 
of Dorey et al. The generalization to the system of Alday et al. is straightforward. Section 
5 is devoted to the analysis of the resulting equations. With strong support from numerics 
we argue that the transition to the mixed regime is associated with the simplest kind of 
excitations in the auxiliary integrable system. Once this is established, the precise value of 
E<i can be calculated analytically. A few technical details are contained in two appendices. 

2 Multi-Regge kinematics of six-gluon amplitudes 

We consider the scattering of six gluons with momenta pj,j = 1 , . . . , 6 satisfying the on- 
shell condition p 2 = 0. From these momenta we can form Lorentz invariant combinations 
of the form 

Throughout this note we shall use the metric r\ =diag(—, +, +, +). Furthermore, we extend 
the range of the index j on pj to the integers with the periodicity condition pj = Pj+Q. 
Due to momentum conservation the quantities xfj are symmetric in their indices, i.e. the 
obey 

— ^ Six point amplitudes of a 4-dimensional gauge theory depend on eight 
kinematical invariants. Hence, there exist many relations between the xfj we introduced 
above. 

In the literature, six-point amplitudes are usually parametrized in terms of the following 
nine invariants, as shown schematically in figure 1, 



S — ~ X \% ! 


Si = 


-T 2 
x 24 J 


S2 = 


2 2 
~ x 35 1 s 3 — — ^46 > 






tl = 


r 2 
x 13 > 


t 2 = 


r 2 

X 15 ! 


(2.2) 


S345 = 


r 2 
x 25 J 


S234 = 


r 2 

x 14 J 


2 

•S456 — -^36 • 





Variables denoted by s or t are referred to as energies and momentum transfers, resp. In 
the Regge limit, energies are large (and can be positive or negative), whereas momentum 
transfer variables are negative and finite. As an example, a positive energy s implies the 
negative distance x\§. Since we have only eight independent variables, one of the above 
can be eliminated with the help of the Gram determinant, namely det([pj -Pj]i j=l,... 5) = 0. 
Since each of the products [pi ■ pj] may be expressed through the kinematical invariants 
defined in eqs. (2.2), we can use the Gram determinant relation to eliminate one of these 
nine invariants. 
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In the physical regime, all the 's-like' kinematical invariants are positive while the 
't-like' variables are negative, i.e. 

S,Si,S45 6 ,S 3 45 > and ti,S 2 34 <0 . (2.3) 

Finite quantities in a theory with dual conformal invariance can only depend on cross- 
ratios. For a six-point amplitude, there are three independent conformal cross-ratios which 
are defined as 

2 2 2 2 2 2 

_ X 13 X 46 _ ^l g 3 _ x 24 x 51 _ g 1^2 _ x 35 x 62 _ S 2 S 

1L\ 9 9 5 ^2 9 9 3 ^3 9 9 

X 14 X 36 S 234S456 X 25 X 41 5345^234 X 36 X 52 S456S345 

(2.4) 

For later use we have expressed them both in terms of the general kinematical invariants 
Xij and in terms of the nine special invariants that were defined in eq. (2.2). 

In the next section we shall look at the multi-Regge limit for six-point functions which, 
by definition, explores the regime 

s > si,S2,s 3 » —ti, — i 2 , -S234 ~ fixed. (2.5) 

The conformal cross-ratios are functions of the kinematical invariants and one can study 
there behavior in the Regge regime. It turns out that 

U\ ~ U 2 ~ 1 — ^3 — > (2-6) 

vanish in the same order (see Appendix A). Hence, the leading asymptotics of any conformal 
invariant of the momenta pi may be parametrized by the two quantities 

ui = u\/{l — u 3 ) and u 2 = u 2 /{l-us) . (2.7) 

In this simple multi-Regge limit, the six gluon scattering amplitude actually agrees with 
the answer that is predicted by the BDS formula. This is because the limiting values 
, U2,U3)R eg ge — (0,0,1), corresponds to a special case of the collinear limit, namely 
u\ = 0, U3 = 1 — U2 where u 2 — > 0. In the collinear limit, the six-point amplitude reduces 
a 5-point amplitude which is known to be BDS exact. Hence, the Regge limit (2.5) of the 
BDS formula for six gluons needs no corrections and therefore is of no further interest to 
us, at least for the moment. 

As we have outlined in the introduction, however, there exists another Regge limit in 
which the BDS formula requires corrections, i.e. the remainder function does not vanish. 
Reaching this limiting regime requires to pass to a new kinematic region which is defined 
as 

s,s 2 > and s x , s 3 , S345S456 <0, h,t 2 ,s 2 34 <0 . (2.8) 

In the corresponding multi-Regge limit all energy variables are much larger than the mo- 
mentum transfer variables, 

s » s 2 , -si,-s 3 > —tx, -t 2 , -S234 ~ fixed. (2.9) 
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Figure 1. Kinematic configuration for the multi-Regge limit before (on the left) and after (on the 
right) the analytical continuation with the kinematic invariants from eq. (2.2). Momenta of are 
denoted by p, while dual coordinates of the Wilson loop by Xi. 

In order to reach this regime, one either starts from the region where all invariants 
are negative and then continues in s and S2 or, alternatively, one first goes to the physical 
region and then continues all the following energy variables analytically along the following 
simple path 

51(99) = e l,p si, s 3 (ip) = e lip S3 and s 456 (v?) = e^s 456 , s 345 (^) = e lv s 3i5 , 

(2.10) 

where ip E [0, ir] and all the other kinematical invariants are left unaltered. In the process 
of the continuation, the conformal cross ratios u\ and 112 are left invariant while u 3 behaves 
as 

C u : u 3 (ip) = e~ 2i ^u 3 . (2.11) 

At the final point ip = tt all three cross-ratios come back to their original values. Along the 
way, the cross-ratio 113 moves along a full circle. 

Our analysis below starts with some conformal amplitude A{ui) in the physical regime. 
Its precise form will be described in the next section. In a first step we want to continue 
this amplitude analytically along the path (2.11). Since A(ui) is not single valued in the 
complex u— space, the resulting amplitude A c (u{) differs from the amplitude we started 
with. Then we perform the Regge limit. This leaves us with a function a{u\,v,2) that 
describes the leading Regge asymptotics of A c . It is this quantity we are actually after. 

The amplitude we shall subject to this analysis is taken from the recent work of Alday 
et al. By construction, it is initially defined in the Euclidean region in which all the 
kinematical invariants (2.2) are negative. In order to go to the physical region (2.3) one 
needs to continue s,Si,s 3 4^ and S456 analytically while keeping ti and S234 fixed. This 
continuation may be performed such that all three conformal cross-ratios ui remain fixed. 
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Hence, conformal amplitudes do not require analytical continuation in order to pass from 
the Euclidean to the physical region. 



3 Six gluon amplitudes at strong coupling 

In the first part of this subsection, we shall briefly recall the prescription of Alday et al. 
for six gluon amplitudes at strong 't Hooft coupling. The second subsection contains an 
analysis of the Regge limit in the physical regime. We shall demonstrate that the Regge 
limit of the six gluon amplitude is BDS exact, in agreement with general arguments. 

3.1 Review of results from Alday et al. 

According to the work of Alday et al. [28], the strong coupling limit of a six gluon amplitude 
in Af = 4 supersymmetric Yang-Mills theory is obtained through exponentiation of the 
quantity 

4(6) . 4( 6 ) i 4( 6 ) 1 4( 6 ) 4. 4( 6 ) (r> i\ 

A ~ ^div + ^BDS-like + ^free + ^period ' \ 6 - L > 

at least up to a constant that is independent of the kinematics. The various pieces on the 
right hand side are all known explicitly. In particular the divergent contribution is given by 

4ii = lE l ^(e 2 xl +2 ) , (3.2) 

i 

where e is the cut-off. The BDS-like amplitude is also relatively simple. In the case of six 
gluon scattering, it is given by 

^BDS-like = -l^[ l °S 2 xl i+2 + ^(-l) k+1 logxl i+2 logx^ +2k+li+2k+3 ) . (3.3) 

i=l V k=0 / 

For us, it is most relevant to have an expression for the difference between the usual BDS 
and the BDS-like amplitudes. This difference is given by 



a (6) = 4 6 i s -4 6 kii ke = EQiogV + ^u-.,)) • 



(3.4) 



The conformal cross-ratios Ui that appear in this expression were defined in the previous 
section. 

The two remaining contributions to the total amplitude cannot be written directly in 
terms of the three cross-ratios Instead, they involve a new set of parameters m,C,(f) 
that are related to the Uj through some functional relations to be spelled out below. Given 
m, C, (j) Alday et al. instruct us to solve the following system of non-linear integral equations 
for three functions Yi = Yi(6),i = 1,2,3, 

/oo 
dO'Ki{0- 0')log(l + Y 2 {9')) 
-00 

/oo 
d9'K 2 (9 - 9') log((l + + Y 3 (9'))) , (3.5) 

-00 
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/oo 
dff 'K 2 (9 - 9 ')log(l + Y 2 {9')) 
-oo 

/oo 
d^ 1 (0-Olog((l + il(O)(l + >3(^))) • (3-6) 
-oo 

These formulas can be used as long as | Im(# — i(j>)\ < ir/2. For values of 9 outside this 
strip, one can employ certain recurrence relation to determine Yf, details are given in the 
Appendix B. The equations (3.5) and (3.6) involve two integral kernels K\ and K 2 that 
are given by 

Suppose for the moment that we can solve these equations somehow. The resulting 
functions Yi{9) clearly depend on the parameters m,C and <j> in some complicated way. It 
is now possible to express m, C, (f> through the cross ratios Ui by solving the following three 
equations 

X2 ™>2 / i 

1 / - -L 



U\ — —r 5- — 11 + 



AA V Y 2 (6 = -in/4) 

"■ = ^ = ( 1+ nra)"- (3 - 8) 

x 35 x 26 ( ~\ 1 1 



3 4 6 ^5 V + r 2 (0 = -3m/4) 

The last equation contains the function Y 2 outside the strip \Im{6 — i(f>)\ < ir/2. With the 
recurrence relations described in the Appendix B we rewrite this equation in the following 
form 



(1 - u 2 )u 2 

U3 - 



(1 + Yi)u 2 + (« 2 + Y 1 )Y 3 



(1 - ui)ui 



(1 + Yi) Ul + (ui + Yi)y 3 



(3.9) 



Once m = m(ui),C = C{uj) and <f) = cp(ui) are known, we can construct the two missing 
terms of the total amplitude. The so-called free energy is given by 

A hl = / ^mcosheiog[(l + Y 1 {9 + i(j)))(l + Y 3 (e + i(l)))(l + Y2{9 + i(l)))^j . (3.10) 

The final contribution to the amplitude depends on the cross ratios only through the pa- 
rameter m = m(ui), 

4eL = \™ 2 - (3-11) 
This concludes our description of the input from the work of Alday et al. 

3.2 Expansion for large masses 

It is instructive to first study the Regge limit in the physical region and to verify that the 
remainder function does indeed vanish (up to a constant). Most of the formulas in this 
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subsection have appeared in the literature before, see [32, 33]. We repeat them here mostly 
for pedagogical reasons. 

In the region (2.3) one can show that the multi-Regge limit is mapped onto large values 
of m, while C remains constant and <f> goes to zero. To see this we start from eqs. (3.8), 
(3.9) and observe that our Regge limit (2.6) forces the values Y 2 (±iTr / 4) , Yi(0), Y 3 (0) to 
approach zero at the same rate, i.e. the right hand sides of eqs. (3.5) and (3.6) must become 
large and negative. As long as \4>\ < 7r/4, the parameter m is forced to become large and 
the first terms on the right hand side dominate, i.e. 

log Y 2 (60 ~ -W2cosh((9 - uf>) , (3.12) 
logF 2±1 (0) « -m cosh(# - i<f>) ± C . (3.13) 

As we shall demonstrate further below, these approximations are justified since the integrals 
on the right hand side of eqs. (3.5) and (3.6) remain finite. This will no longer be that case 
after analytic continuation into the kinematic region (2.8). Taking the difference of the two 
equations (3.6) we obtain 

Y 3 (0)/Y 1 (0) = e 2C . (3.14) 

As we argued above, the Regge limit requires both Y\(0) and Y^(0) to approach zero at the 
same rate so that their ratio stays fixed. Hence, in our Regge limit C is a finite constant. 
Finally, the ratio of u\ and U2, which should also be nonzero and finite in the Regge limit, 
is given by 

ui/u 2 ~ Y 2 (9 = -m/A)/Y 2 {9 = ivr/4) + . . . . 

Here we have omitted terms of lower order in the Regge limit. Using eq. (3.12) we conclude 
that msmcj) has to be finite, i.e. that (j) tends to zero. 

For the following analysis of the Regge limit it is convenient to introduce the parametriza- 

tion 

£ = e -mcos0 ^ w = e msin0 > Q = cosh(C) _ (3 15 ) 

We shall expand all equations from the previous subsection in the limit e — > while keeping 
both w and c fixed. From eqs. (3.8) we find that 

u\ = we + O (e 2 ) , 
u 2 = w~ 1 e + 0(e 2 ) , 

u 3 = l-(w + 2c + vj- l )e + O {e 2 ) . (3.16) 

One easily verifies that the dependence of the conformal cross ratios on e agrees with their 
behavior in the Regge limit, i.e. u\,u 2 and 1 — u 3 vanish in the same order, when e is 
sent to zero. The leading contributions to all conformal invariants in the Regge limit are 
parametrized by the two functions 

itl w 2 _ u 2 1 

1 1 — ^3 w 2 + 2cw + 1 2 1 — ^3 w 2 + 2cw + 1 
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We can invert these two equations to determine w and c as functions of u\ and U2, 

« W— , c(iti) « . (3.17) 

V «2 2y/U\U2 
Our aim now is to evaluate the dependence of the so-called remainder function R, 

- flWfc) = AW - A£ - = A^ + Ag) iod - A® + constant , (3.18) 

in the Regge limit as a function of U\ and U2- The function A was spelled out in eq. (3.4). 
In our limit we find that it is given by 



_ A(6) _ A (6) _ A6) 

~ ^BDS-like BDS 



i=l 



= -\ \og\w) - \ log 2 ( £ ) - g + O (e) . (3.19) 

The divergent contribution is actually cancelled exactly when we subtract the period con- 
tribution to the remainder function 

^lod = \™ 2 = ^log 2 H + ^log 2 (,) . (3.20) 

It now remains to compute the contribution from the free energy to leading order in the 
Regge limit. When we send m to infinity, the functions Yi become small all along the real 
axis so that we can approximate log(l + Yi(8)) ~ Yi{6). In this limit, the free energy then 
takes the form 

A t!L ~ m j y cosh eY i( e + i( t>)+ m j Y cosh 0Y 3 (e + i(t>) 

/df) 
— cosh 9 Y 2 (9 + i(f>) . (3.21) 
27T 

All three integrals can be evaluated easily with the help of the following integral formula 

/OO ]Q 1 
— cosh ee- amcoshd = -K x (am) , (3.22) 

such that 

A hl ~ — coBh(C)^i(m) + — lfi(V2m) . (3.23) 



7T 7T 



In the limit of large m we can use the asymptotic form of xK\(x) ~ w nx/2 exp(— x). Since 
the parameter m depends on e through 

/i 2 , i 2 , llog 2 W 

m = (log £ + log w I loge h . . . 

v ' 2 loge 
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we conclude that the free energy vanishes when e is sent to zero. Indeed, once we insert 
the leading dependence of m ~ — loge in e into the formula (3.23) for the free energy, we 
obtain 

4 6 i = /f ce^g-e + ... . (3.24) 

where . . . represent lower order terms. Combining this result with our previous equations 
(3.19), (3.20) and (3.24) for the Regge limit of the amplitude A and the periodic contribution 
we obtain, 

-R {6) (ui) » -fj + constant + 0(e) . (3.25) 

As we had anticipated, the leading Regge asymptotics of the remainder function in the 
physical regime (2.3) is a constant. 

4 Excited states in thermodynamic Bethe Ansatz 

In this section we shall perform the analytic continuation of the kinematic variables along 
the path (2.10), and we show that excited states of the system start to contribute. To explain 
the basic ideas, we begin the presentation of a toy model [30, 31]. The corresponding analysis 
of the system (3.5) and (3.6) is more involved and is performed in the second subsection. 

4.1 Analytical continuation and excited state TBA 

As a warmup example we want to study a rather simple non-linear integral equation that 
depends on a single complex parameter m. We are interested in the behavior of the energy 
Fo(m) as it is continued analytically in the parameter m. Our toy system involves a single 
function Y(8) which is determined by 

1 f°° 

- log Y{9) = m cosh 9 / K{9 - 9') log (l + Y(9')) ■ (4.1) 

We shall think of the kernel function K(9) as being determined through an S-matrix, 

K{9) = -i^ log 5(0) . (4.2) 

If the S-matrix is unitary, i.e. S(9)S~ 1 (—9) = 1, then the function log S(9) is antisymmetric 
and so K{9) is invariant under 9 — > —9. Let us furthermore assume that the ground state 
free energy is given by 

o poo 

F (m) = — d9 m cosh 9 log (1 + Y(6)). (4.3) 

In order to understand how Fa(m) depends on the mass, it is quite instructive to study 
a trivial system with constant S-matrix. In this case, the equation (4.1) determines the 
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solution logY(#) = —m cosh 9 and we can evaluate the free energy of the ground states 
exactly, 



Fo(m) 



7T 2 . 


roc 

/ dfl 

' —oo 


M 


3m 2 


2 


2vr 2 


6 


OO y- 



d0 m cosh 9 log(l + e" 



-mcosh$\ 



7T 



fc=l 



A: 



log m + - + log 7T - 7s 



+ ^ ( ^m 2 + (2fc - 1) 2 tt 2 - (2fc - I] 



m 



7T 



fe=l 



2(2/c - l)vr 



(4.4) 



with the Euler-Mascheroni constant — lini s ->oo(I]m = i m_1 — 1°S S ) = 0.577215.... The 
branch cut of the logarithm is chosen to lie along the negative imaginary axis. For the 
integration we have used formula (8.526) from [34]. Then we expressed the infinite sum 
over Bessel functions K\ through simpler functions with the help of 



oo 1 oo 

^K (A;x) = -( 7E + log| F ) + 7r^ 



k=l 



k=l 



y/x 2 + (2k 



1 

2kn 



(4.5) 



In addition, we used the formula [x-RTi(x)]' = — xKq(x) for the derivative of xK\(x). The 
free energy Fo(m) possesses a series of square root singularities at the points x^ = i(2k—l)n. 
Upon analytical continuation along a curve C that encircles the first N of these singularities 
in the upper half of the complex plane, we obtain 



F c (m) = F (m) 



12^ 



7T 



fc=l 



yjm 2 + (2k 



1) 2 TT 2 



(4.6) 



We shall reproduce this simple answer down below through an analysis that does not require 
the S-matrix to be trivial. 

When the system possesses a non-trivial S-matrix, the free energy can no longer be 
evaluated exactly. Nevertheless, we can still analyse how the the free energy Fo(m) changes 
under analytic continuation in the mass parameter m. We begin by re-writing our non-linear 
integral equation (4.1) in the form 



io g y(0) 



m cosh 9 

2vr 



m cosh 9 H 

2vr 



d9'(-i^- Q log S(9 - 9')) log(l + Y(9')) 



d9'i log S(9-9') 



Y'(9') 
1 + Y(9') 



(4.7) 



Note that the integrand has poles at those points #o at which the solution Y(9) assumes 
the value Y(9q) = —1. Since the solutions Y(9) depend on the parameter m, the solutions 
of Y(9q) = — 1 depend on m, i.e. they move through the complex plane as we vary the 
parameter m. Whenever one of these solutions approaches the integration contour, we 
deform it so that solutions never cross. We do that until we get back to the original value 
of m. At this point we now determine a solution Y c (9) through an equation of the same 



- 12 - 



form as eq. (4.1) but along some oddly shaped contour that depends on C. We want to 
move the contour back to the real axis. Along the way we pick up contributions from the 
solutions of the equations 

Y c {6i) = -1 . (4.8) 

Once the integration contour is back on the real axis, our non-linear integral equation takes 
the form 

-lo g y c (0) = mcosh0- — / dO'K{9 -6')\og(l + Y c (6')) - ^(-l)"Mog £(0 - 0;) , 

2?r J-oo i=1 

(4.9) 

with integers rii depending on the contour. The two equations (4.8) and (4.9) determine 
the function Y c along with the parameters 0j. Similarly, from the following formula for the 
free energy 

o [-co o rco Y'(9) 
F (m) = 77 / d0m cosh 01og(l + r(0)) = k / d0msinh0 ^-TTTT ■ ( 4 -!0) 

^ J-oo 7T 1 + Y(6) 

we infer that 

F c (m) = -s / cZ0mcosh01og(l + Y c (0)) V (-l) n H sinh 4 . (4.11) 

In order to compare with our previous result (4.6) for the trivial system, we set S = 1 
so that — logY c = — logy = mcosh0. Suppose that m is continued along a curve in the 
upper half of the complex m-plane that encircles the points m = iir(2k—l) for k = 1, . . . , N. 
Then there are 2N solutions of the equation Y(0) = — 1 that cross the real line. Since these 
are located at ±0^ with mcosh0fc = ±i(2/c — l)ir, formula (4.11) becomes 

12m N 12 - 
F c (m) = F (m) Visinh0 fc = F {m) > ^f{2k - 1) 2 tt 2 + m 2 . 

k=l k = l 

Here, we expressed sinh0fc through cosh0fe. The sum extends over ./V pairs of solutions ±0&. 
Note that both solutions within each pair contribute the same amount. Therefore, we have 
simply multiplied the sum over iV terms by an overall factor of two. Our result is identical 
to our previous formula for Fq in a trivial system. 

We conclude our analysis of the toy model with a few important comments. The 
simplest non-trivial case of a continuation occurs when a single pair ±0o of solutions to 
Y(9q) = — 1 crosses the real axis. In that case, our equations for Y c ,Fq and 9q become 

-logY c (0) = m cosh + log S Jf - 6 o °\ - -L r d6'K{e - 9') log(l + Y c {6')) (4.12) 

o(0 + o ) 2tt 7„ 00 

F$(m) = isinh0 o + ^ / d6 m cosh log(l + Y c (6)) . (4.13) 

The position of the 6q can be calculated making use of 

-lo g y c (0 o ) = i(2Jfe-l)7r = mcosh0 o - logS(20 o ) +log5(0) 



i r°° 

— / d0 / i^(0 o -0 / )log(l + y C (0 / )). (4.14) 
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The first equality is obtained from eq. (4.8) by taking the logarithm. The second equality 
just repeats the functional equation for Y c . According to the previous equation, the right 
hand side has to stay finite when we send the mass parameter m to infinity. Since the first 
term on the right hand side diverges with m, this divergence must be cancelled by one of 
the other terms. 2 In the limit m — > oo, the integral is a small correction and <S(0) does not 
depend on m at all. Hence, the term log<S(#o) has to diverge with m, i.e. the parameter 
20q has to approach one of the poles 6s of the bare S-matrix. As we go to finite values of 
the mass parameter, this value of 9q will receive corrections, 

6 = s /2 + 5 . (4.15) 

In principle, we can compute such corrections perturbatively in the mass parameter m. It 
will turn out later that the case of two poles crossing the real line is particularly relevant 
to the study of the Regge limit. 

Before we return to the study of the equations (3.5) and (3.6), let us make one comment 
on the path of the analytic continuation. In this subsection we assumed the curve C to 
be closed in the complex m-plane, i.e. in the parameter space of the non-linear integral 
equation. This is important for the relation with excited states where m plays the role of 
a physical mass parameter. Our derivations, on the other hand, remain perfectly valid for 
arbitrary curves C that begin at m and end at a possibly different point m! of the m-plane. 
All we have to do is to replace m by m! in our equation (4.9) and all equations after eq. 
(4.11). This extension shall play an important role in the subsequent analysis. 

4.2 Excited state TBA for six gluon amplitudes 

It is not difficult to extend the analysis of the previous subsection to the system that 
describes scattering amplitudes of six gluons. In this case, we need to solve the equations 
(3.5) and (3.6) for three functions Yj, i = 1, 2, 3. The form of these equations is very similar 
to eq. (4.1), with one notable difference. In the context of six-gluon amplitudes the physical 
parameters are the cross ratios U{ rather than the parameters m, C, (j). Since these two sets 
of variables are related by some non-linear equations (3.8), an open curve in the space of 
parameters m, C and eft may be mapped to a closed curve in the cross ratios. Consequently, 
we shall allow our curve C to end in a point m',C' and <fi' that may differ from starting 
point m, C and <fi. Dealing with the requirement u[ = U{ is differed to the next section. 

When we analyzed our toy model is was crucial that the kernel function K could be 
obtained from an S-matrix through equation (4.2). In complete analogy, we need to find 
unitary, crossing symmetric S-matrices 

^ log Sj(9) = -27^(0), (4.16) 
for the two kernels K\ and K<i that were defined in equation (3.7). The solution is given by 

Sl(9) ~ 'TT^' S2(9) ~ 2*sinh(0) + V2 • (417) 



2 Here we assume that the solution of eq. (4.8) does not coincide with the zeroes of cosh#o- This holds 
true for the cases to be considered in the next section. 
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Up to a common sign, the solution is uniquely fixed by the unitarity condition 

S(9)S(-9) = S(9)S*(9*) = 1. (4.18) 

One may also verify that Si are crossing symmetric, i.e. that they obey 

Sj(6) = S ai (9) = S ab (iir-6) = Sj(iir-0), (4.19) 

where S 2 (0) = S 2 (0) and S\{9) = —S\{9). Once we have rewritten the functional equations 
in (3.5) and (3.6) in terms of Si, we can perform a partial integration and analyze the 
behavior of solutions Yi upon analytic continuation of the parameters m, C and (j). The 
results are 

/oo 
d9' K^9 - 9')log{l + Y 2 C [9')) 
-oo 

Ni 

d9'K 2 {9 - 9') log((l + if (0'))(1 + Yi(0'))) + ]T(-ir- log S 2 (9 - 0^) 

3=1 

-2 ^T(-l) n ^ log S 1 (9 - 9 2J ) - £(-l) n3 - log S 2 (9 - 9 3J ) , (4.20) 

3=1 3=1 

and similarly 

/oo 
d9'K 2 (9-9')log(l + Y 2 c (9')) 
-oo 

/oo N l 
d9'K x {9 - 9') log((l + Y?{e')){l + Y 3 C (9'))) - ]T(-1)^ log S^O - ltj ) 

3=1 

N 2 N s 

- ^(-lp-' logS 2 (9 - e 2J ) - log S 1 (9 - 6 3tj ) , (4.21) 

3=1 3=1 

with integers N s and n s j depending on the path of the analytical continuation. The pa- 
rameters 9 s j are obtained by solving the equations 

Y S C (9 S , 3 ) = -1 . (4.22) 

As before, we can now evaluate the formulas (4.20) and (4.21) at the points 9 s j and insert 
the conditions 

-logY s c (9 sJ ) = i(2fc sJ -l)7r (4.23) 

to obtain equations for the positions 6 s j. Since the resulting expressions are rather lengthy, 
we shall refrain from spelling them out. Similarly to the toy model considered in the 
previous subsection, the position to the nearest poles in the limit of large m are related 
to poles of the S-matrices Sj. Eqs. (4.23) and (4.21), for example, require that for every i 
at least one of the differences 9\^ — 9 2 ±\j or (?x,i — $2,j must approach a pole of S± or S 2 , 
respectively. 
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The role of the free energy Fq is played by the functional that was defined in eq. 

(3.10) . During the analytical continuation, this quantity changes according to 

Af£ = j ^m' cosh 9 log [(1 + Yf(0 + i(j>')){l + Y 3 C (6 + i<f/))(l + Y 2 C {6 + itf))^ 

Ni N3 N 2 

-m'i^(-l)™ 1 ^ sinh^ij - m'i^{ 1 j "•'••< sinh# 3ij - yftm'i ^{-lf 2 ^ sinh# 2ij , (4.24) 
j=i j=i j=i 

where variables S stand for S j = S j — i<jJ '. The possible s j that contribute terms 
to the second line are solutions of eq. (4.23) that have passed the line Im# = —<p upon 
analytic continuation. For nonzero cf>, the number of such solutions and the direction in 
which they cross the line may differ from the corresponding numbers and directions for 
the real line Im# = 0. This is why we had to use new integers N s and h s j instead of the 
corresponding parameters N s and n s j that appeared in the integral equations (4.20) and 
(4.21). With these comments we conclude our general discussion of analytic continuations 
in the parameters for six gluon scattering. We shall now return to the Regge limit. 

5 Regge limit of six gluon amplitude 

Now we are prepared to combine the results of the previous sections in order to complete 
our main goal, namely to determine the Regge limit of the six gluon amplitude in the 
regime (2.8). As explained in section 2, the first step involves an analytic continuation in 
u-space. The transformations (3.8) describe m, C,(f) as a function of the conformal cross 
ratios Uj and hence they determine how the parameters of the integral equations (3.5) and 
(3.6) change under analytic continuations along C u . Some important features of the curve 
in the space of parameters m, C, (j> are discussed in the first subsection. Numerical results 
show that very few solutions of the equations (4.22) actually cross the real line while we 
continue from the physical to the mixed regime. This enables us to determine an explicit 
formula for the Regge limit of the six gluon amplitude in the regime (2.8). 

5.1 The path of the analytical continuation 

Our first aim is to describe the path of the analytical continuation. As we explained in 
section 2, we would like to continue analytically along the path C u that was defined in eq. 

(2.11) with u\ and 112 kept fixed. From the first two relations in eq. (3.8) we infer that 

Y 2 - d ^ y 2 Hvr/4) = Y+ d 4 f = . (5.1) 

1 — U\ 1 — U2 

Note that the left hand side depends on the parameters m,C,(f> which enter the nonlinear 
integral equations for Yi. Since we want to keep the right hand side fixed, eqs. (5.1) provide 
two constraints on the form of the curve in the space of parameters m, C, (j). 

We still need to explore one more relation that involves the cross ratio U3. Let us recall 
that U3 has been expressed in eq. (3.9) through the values of the functions Y{ in the strip 
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Figure 2. Upon analytic continuation along C, the parameter m is shifted from its initial value. 
The plot shows numerical results for m = 10, coshC = 3/5 and = at ip = 0. Plots for different 
values of the parameter C(tp = 0) differ by terms of order exp(— m). 



\Im(9 — i(p)\ < 7r/2. We can actually turn eq. (3.9) into a formula for C = C(ui) once we 
notice that 

*3(0)/n(0) = WVIMO) = e 2C . (5.2) 

This simple relation follows by subtracting the non-linear integral equations (3.6) for Y\ 
and Y3 from each other. We employ eq. (5.2) to eliminate both Y3 and Y\ from the equation 
for M3 and arrive at, 

cosh 2 C = (l-"i-"2-u 3 ) 2 , (5 . 3) 

4uiU 2 M3 

During our continuation we keep u\ and U2 fixed while U3 is moved along a full circle 
until it comes back to its original position. Because the left hand side involves the square 
of c = coshC, the quantity coshC has changed its sign once we reach the end-point of the 
curve C u . Note that eq. (5.3) determines the parameter C explicitly as a function of the 
cross ratios. The variation of the real and imaginary part of C along the path C u is plotted 
in fig. 3. All the curves start in the lower half on the left hand side and proceed to the 
central line at ImC = ir/2 and return to the left and into the upper half of the figure. At 
(p = 7r/2, the real part of C is of order m. 

The dependence of m and cj) on (p is a little more difficult to obtain. It requires inverting 
the eqs. (5.1) which contain the function Y2. The latter must be determined numerically 
by solving the non-linear integral equation for a large number of values m, (ft and C. Given 
Ui we can restrict the values of C to 

ctfcfc.) = "-";-"~ 2 " >2 , 
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Figure 3. The plot shows the dependence of C on (p. The various curves correspond to the values 
m = 10, = and coshC = c at the starting point tp = 0. Here, we plot C in its complex plane 
along the path given by (5.3) while 4>{<p) = 0. 

where cp runs from (p = to cp = ir. We have performed these numerical studies for special 
cases in which u\ = U2- The condition amounts to setting = at p = 0. It is not difficult 
to see that = will remain true as we move along the contour C u . Hence, it remains to 
determine m(ip) such that l2(=t* 7r /4), and hence the cross ratios u\ = U2, remain invariant 
as me change <p. An example is shown in fig. 2. The shape of the curve has very little 
dependence on the parameter C. We observe that the parameter m becomes complex. At 
ip = tv, its imaginary part is Imm(7r) = — tt. The real part also gets shifted by an amount 
Am = Rem(7r) — Rem(0) ~ 1.8. We shall reproduce these shifts through an analytical 
argument in the next section. Some more comments on the numerical studies can be found 
in appendix C. 

5.2 Evaluation of the Regge limit 

In the previous section we have investigated the image of curve C u in the space of parameters 
m, (j) and C. For = 0, the corresponding values of m(ip) and C{(p) are depicted in figs. 2 
and 3. These numerical data provide us with a good qualitative understanding of the curve 
along which we would like to continue the parameters of the non-linear integral equations 
(3.5) and (3.6). As we discussed in section 4, we need to watch the solutions of Yi{0) = — 1 
while continuing m, <p and C. Our numerical calculations show that, for large values of 
m, very few poles get close to the integration contour. In the case of I3, the numerical 
results are displayed in fig. 4. The plot was produced for a continuation along the paths 
described by figs. 2 and 3. We see a single pair of poles that crosses the real line. The 
symmetry under reflection — )■ —9 is a direct consequence of our special choice = 0. 
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Analogous results for the solutions of Y 2 = —1 are displayed in fig. 5. In this case, one pair 
of solutions approaches the real line while all others stay away from it. We did not display 
the corresponding plot for Y\ since it looks very similar to that for Y3, except for an overall 
shift by ztiir. 

With these results in mind we can actually spell out a system of equations that should 
determine the solutions Yi at the end of the analytic continuation, i.e. in the mixed regime, 

/oo 
d9'K 2 (0 - 6') log((l + Y((9'))(l + Y^(9'))) 
-oo 

/oo 
d9'K x (9 - 9') log(l + r 2 '(0')) + log S 2 {9 - 9 X ) - log S 2 {9 - 9 2 ) , (5.4) 
-00 

/oo 
d9'K x {9 - 9') log((l + Y[{9')){1 + Kj'(0'))) 
-00 

/oo 
d9'K 2 {9-9')\og{l + Y 2 \9')) +lo g Si(0-0i)-logSi(0-0 2 ) ■ (5.5) 
-00 

One might have expected to see additional terms from the solutions of Y 2 = — 1 on the right 
hand side. As one may easily see, these take the form 

lim(21og(Si(0-a))-21og(Si(0 + a))) = lim (4i sech(0)a) = 0, (5.6) 

a— >0 a->0 

and 

lim(21og(S 2 (0-a))-21og(S 2 (0 + a))) = lim ( 8iV2 cosh(0) sech(20)a) = 0. (5.7) 

a— 5>0 a->0 V / 

Hence, we there are no additional terms on the right hand side of eqs. (5.4) and (5.5). In the 
above equations and throughout the entire discussion below we place a / on all quantities 
that have been continued to the endpoint of the curve C = C u that we described in the 
previous subsection. The parameters 9\ and 2 appearing in the arguments of the two 
S-matrices (4.17) must solve one of the equations 

Y{{e v ) = -1, Yi{9 u ) = -1 (5.8) 

depending on value of C. In addition we shall use the convention that, after crossing the 
contour, 9\ has positive imaginary part while 2 has negative imaginary part. For = 
we have 9\ = —9 2 . 

In the leading approximation for large values of the mass parameter m, the non-linear 
integral equations simplify, 

logY 2 '(0) = -m'V2cosh(0-# / )-log5 2 (0-0i)+log5 2 (0-0 2 ), 

logY 2 ' ±1 (0) = -m / cosh(0-^')±C"-log5i(0-0i) + log5i(0-0 2 ) . (5.9) 

As we explained earlier, the solutions to eq. (5.8) are determined by the poles of the S- 
matrix, at least to leading order in large m. Since S 2 has poles at 9$ with expi9s = i, the 
solutions 9 V ~ 0s/2 are well approximated by 

01 = -02 = i\ + ••• , (5.10) 
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Figure 4. Positions of some central solutions to the equation Y^{ff) = —1 along the path depicted 
in figs. 2, 3. Light dots correspond to positions in the first half of the path before any solution has 
crossed the real axis while the dark ones are associated with the second part. 
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Figure 5. Positions of some central solutions to the equation ¥2(6) = —1 along the path depicted 
in figs. 2, 3. We change from light to dark dots at the point where the first solutions of Y3 = — 1 
cross the real axis. 
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with corrections of order exp(— m). We can insert the formulas for Y- with 9\ = — 62 = i7r/4 
into the equations (5.1) and (5.3) that determine the conformal cross ratios at the end- 
point of the analytic continuation. A short computation gives 

u [ = w ' 7£ > + o ((e') 2 ) , 

u > 2 = (w'y^e' + O^e') 2 ) , (5.11) 
u > 3 = 1 - (y/ - 2c' + (w'y 1 ) >ye' + O {(e') 2 ) , 

with 7 = — (3 + 2\/2). Here, the parameters w',c' and e' are determined by our formulas 
(3.15) only that we have to replace m, <j>, C by m' , <fi' , C . At the end of the analytic contin- 
uation we want to return to the original values of the cross rations, i.e. we need to impose 
Ui = u^. This leads to 

e' = 7" 1 e + 0(e 2 ) , w' = w + O (e) , c = -c + 0{e), (5.12) 

Obviously, there exists a second solution of the conditions m = u^. This solution, however, 
violates the condition c = —d that we have argued for in the previous subsection. 

It is instructive to re-express our relations (5.12) through the parameters m,C and <p. 
We have seen already that coshC = — coshC. For the other two parameters we obtain 

m! = ^m 2 + s 2 + 2ms h cos <j) + O (e" m ) , (5.13) 

m'cos(f)' = Sh + m cos <f> + O (e" m ) , (5-14) 

where Sh = log 7 = — iir + log (3 + 2v / 2). Let us recall that these relations between the 
value of m, C and <j> in the initial and final point of our curve C were derived under the 
assumption of single pole crossing. The numeric results we displayed in the previous section 
were obtained for <fi = 0. In this case, the assumption of single pole crossing gives <j)' = 
and 

m! = m- ivr + log(3 + 2\/2)) + O (e~ m ) . (5.15) 

This coincides with the shift we found in our numerical studies, see fig. 2, and provides 
indirect evidence for the assumption made in the analytic treatment. 

Let us now turn to the main goal, namely to the calculation of the remainder function 

_#(6) = 4(6) + A ' p (6). od _ A /(6) + constant) (5 . 16) 

where the individual terms are not evaluated after analytic continuation along the curve 
C. Let us begin with the function A' that was spelled out in eq. (3.4). Upon continuation 
along the curve (2.11) we find 

A /(6) = A (6) _ ^ (logn3 _ log (l _ n3 )) _ Hi 

2 

= ^irloge + \iTr\og(2c + w + W" 1 ) -y + A (6) + O (e) . (5.17) 
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In passing from the first to the second line we have inserted the asymptotic form (5.11) of 
the cross ratios u' i = U{. The period contribution is even easier to find, 

A'£L = \M 2 = \(mr-\{mf + Af enod 

= - l - log 7 log(e) + \ log 2 7 + ^ 6 e li 0d + O (e) 

= ivrlog(e)-ilog| 7 | log(e) + i log 2 7 + period + ■ ( 5 - 18 ) 

Finally, we also have to determine the contribution from the free energy. Specializing our 
general result (4.24) to the present case we find that 

42 = / ^m' cosh log [(1 + Y[(9 + i<t>')){\ + Yi(9 + icj>'))(l + Yi(0 + 

-m'i sinh(0i - i<j>') + m'i sinh(6> 2 - i<j>') + V2m' sinh(|0'|) . (5.19) 

Note that the second line contains a contribution from one solution of the equation Y 2 = — 1. 
As we explained at the end of section 4.2, the contributions to the free energy arise from 
solutions 6* that cross the line Im(9* = <j). As one can read off from fig. 5, one of the 
solutions of Y2 = —1 does cross this line, at least when (j) ^ 0. 

In the limit of large m, the integral in eq. (5.19) gives a vanishing contribution as one 
may see by repeating the analysis we have carried out in great detail in section 3.2. So the 
only terms that are relevant in the Regge limit come from the second line of the previous 
equation. This gives 

AfJ e = V2m' cos(^) + V2m' sinh(|^|) + O (e) = v^log e' + y/2\ log w'\ + O (e) 

= -\/21oge + \/21og7 + 42e + logio| + O (e) . (5.20) 

In passing to lower line, we have inserted eq. (5.14) and also added the term Af ree that was 
evaluated in section 3.2. Since the free energy before analytic continuation vanishes in the 
limit of e — > 0, the additional term is actually zero. We only put it back in to show that 
the leading term should be regarded as a difference of energies. Using the final result (3.25) 
of our computation in section 3.2, we are now prepared to determine the Regge limit of the 
remainder function in the regime (2.8), 

_ D'(6) _ 4/(6) 4/(6) _ A /(6) 
U — ^free + ^period ^ 

= -e 2 log e - \i log (2c + w + w' 1 ) + y/2\ log w\ + constant' + O (e) . (5.21) 

Here, we have moved all the constant terms into a redefinition of the 'constant' and defined 
the constant e 2 by 

e 2 = (y/2 + \ log(3 + 2v^)) . (5.22) 

In a final step, we want to express the variables w, c and e through the conformal cross 
ratios. Using eqs. (3.17) along with 



V%^2( 1 ~ u 3) , (5-23) 



- 22 - 



we find that 



R' = -e 2 log((l - n 3 )v / tiitt2) + ^inlog y 'u\Ui + V2\ log(ui/n 2 ) 
-f-const' + O (1 — S3) . 



After exponentiation of this result we obtain 
Amp' ~ (W ) ~ exp 



2vr 



with 



e 2tt ~ e 2 4n 



mtfiifia) ((l-«3)v^ 



exp 



^e 2 
2vr e2 



" "^~T (^div + ^BDS 



R') 



(5.24) 



(5.25) 



U2j \U 2 J \UlJ \Ul 



(5.26) 



up to some irrelevant constant prefactor. This formula is the main result of our analysis. 
It describes the Regge limit of the six-gluon amplitude in the mixed regime and at strong 
coupling. As we have stressed before, our calculation assumes that the analytical contin- 
uation of the amplitude may be performed within the prescription of Alday et al., see our 
discussion in the introduction and the next section. 

Before we conclude, we would like to compare our result with the form of the answer at 
weak coupling. The relevant expression is known in the leading logarithmic approximation 
[19]. The following formula gives the correction to the BDS expression: 
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(5.28) 



All relevant notations are explained in [19]. In particular, the % are the transverse momenta 
of the t{ invariants, the parameter a is given by a = A/87T 2 and A is the 't Hooft coupling 
A = g 2 N c . The exponent w of s 2 is given by 



/ Tl Tl \ 

u(u,n) = Aan^{l)-^{l + iv+-)-i>{l + w--)j 



(5.29) 



In the weak coupling expression (5.27) the new term, A 2 ^.4, comes as an additive correction. 
There is no indication that this correction should be written as an exponential. From the 
point of view of Regge theory, it seems more natural to sum over the different singularities 
in the complex angular momentum plane, and it is not clear whether exponentiation would 
still be consistent with Regge factorization. However, the leading logarithmic approximation 
alone is not sufficient to exclude the possibility of exponentiation. 
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The function A2-S.4 may be re-expressed in terms of the cross ratios ttj that we have 
built in section 2. In these variables, (5.28) assumes the form 



n=oo C A / ~ \ %v 

™ = \ EH)" T^^-^-l) | coshnC, (5.30) 



where we have used 1 — uj, ~ l/s2 (appendix A) 

<?2P4 




(5.31) 



and 



f Q2PA l-Ui-u 2 -U 3 I-U1-U2-U3 , „ , . 

cos arg = — — « ^=^= = coshu. (5.32) 

V P5?l/ 2 x /«rU2 2^/uiu 2 u 3 

The leading power of S2 belongs to n = 1 and i/ = so that the intercept at weak 
coupling becomes 

w(0,l) = -E 2 = -4(21n2-l). (5.33) 

Consequently, up to logarithmic corrections, the leading asymptotic behavior of the remain- 
der function at weak coupling is given by 

A 2 _> 4 ~ ^(l-na)-^ ' 1 ) coshC. (5.34) 

In this leading logarithmic approximation, it is not possible to determine the scale of the 
energy S2- Put differently, the factor 1 — 113 may be accompanied by some undetermined 
finite function of u\ and U2- 

Let us attempt a first comparison of our strong coupling result eq. (5.26) with the 
weak coupling amplitude in eqs. (5. 27), (5. 28) and eq. (5.34). In both cases, the energy 
dependence of the correction is power-behaved (i.e. Regge-behaved) . We first face the fact 
that, at weak coupling, the correction is additive whereas, at strong coupling, the remainder 
function in the exponent leads to a multiplicative correction to the BDS amplitude. As a 
possible interpretation, we might assume that in eq. (5.27) the correction, which at weak 
coupling is of order A, grows with A and, at strong coupling, dominates over the 1. We are 
then led to confront A2->4 in eqs. (5.30) and (5.34) with eq. (5.26). 

A striking observation is the change in the sign of the exponent E2: The intercept of 
the energy S2 which on the weak coupling side is positive appears to turn negative when 
the coupling is increased. Also, the strong coupling expression exhibits an energy scale 
depending upon the product of u\ and U2- Next let us look at the ratio u\/v,2 which is 
raised to power iv on the weak coupling side. Since at large S2 the leading behavior comes 
from iv = 0, the dependence on U1/U2 drops out of eq. (5.34). In contrast to this, on the 
strong coupling side, the ratio U1/U2 comes with an exponent of order In addition, 

the factor cosh C that multiplies the answer at weak coupling, is absent at strong coupling. 
Finally, in eq. (5.26), we see a phase factor which also depends upon the product of U\ and 
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U2- It should be kept in mind, however, that the full BDS amplitude at strong coupling, 
-^^"'"'^BDS' wnen continued into the mixed region, also contains a phase: A straightforward 
calculation shows that this phase exactly cancels the phase in eq. (5.26). As a result, the 
strong coupling result has no such phase factor. 3 

Let us try a further interpretation and assume that, at strong coupling, the correction 
factor exp(\/XR' /2tt) can also be written in the form of eq. (5.30) as an expansion in the two- 
dimensional conformal group (Mobius group), with some strong coupling analogue 6j(is, n) 
of the function u){y, n). Then, in eq. (5.26), the exponent of u\/u2 could be interpreted as 
coming from a dominant v- value v s that is of the order VX At the same time, the exponent 
of S2 would come from the function u){y, n), evaluated at the value v s . Finally, the absence 
of the factor coshC might result from the fact that, at large u s , ui(v s ,0) dominates over 
uj(v s ,1). In pursuing such an interpretation one would perform a saddle point analysis of 
the ^-integral: At weak coupling, the large parameter is In S2 which leads to the dominance 
of the point v = 0. In contrast, at strong coupling it is y/\ which plays the role of the large 
parameter, leading to a value v<i = 0{yX). In this sense, the form of our strong coupling 
results seems perfectly consistent with the general features of Regge theory. 

6 Conclusions and Outlook 

In this paper we have studied the Regge limit of the planar six-point function at strong 
coupling. In the weak coupling regime this high energy limit, taken in a very special 
kinematic region ('mixed region') where some energies are positive while others are negative, 
is known to be particularly sensitive to the analytic properties of scattering amplitudes. 
Making use of the set of equations for Y-functions, we have been able to perform the analytic 
continuation of the scattering amplitude at strong coupling into this kinematic region and to 
calculate new corrections which are given by excitations of the auxiliary quantum integrable 
system. A comparison of the strong coupling result with the expressions for weak coupling 
shows interesting features but needs further investigations. At the moment, nothing is 
known about how the quantity E2 may be interpolated between weak and strong coupling. 

As we pointed out in the introduction, our calculation of the parameter E2 at strong 
coupling was based on one important assumption, namely that we can simply continue the 
prescription of Alday et al. into the mixed regime. In other words, we have assumed that 
the amplitude in the mixed regime is determined by the exponential of the free energy 
in the auxiliary quantum integrable system. One may certainly wonder whether this as- 
sumption was justified. As we sketched in the introduction, scattering amplitudes at strong 
coupling do possess a geometric origin. According to the original prescription of Alday and 
Maldacena, they are determined by the area of a 2-dimensional surface Sq which stretches 
to the boundary of AdS^. This surface emerges as the unique saddle point that dominates 
the string theoretic path integral in the limit of strong 't Hooft coupling. The kinematic 
invariants of the scattering process are encoded in the shape of the 1-dimensional boundary 

3 In the weak coupling regime, it was the appearance of such a phase factor in the BDS amplitude which 
played a vital role in understanding the deficiency of the BDS formula. The complete cancellation of the 
phases at strong coupling therefore provides support for the correct analytic structure at strong coupling. 
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OSq of Sq in the boundary of AdS$. Our analytic continuation in the space of kinematic 
invariants amounts to a continuation in the boundary conditions for a (path) integral. In 
any saddle point approximation, the contributing saddle points can be very sensitive to 
boundary conditions. In fact, there could even exist points, i.e. values of the kinematic 
invariants in our case, at which several saddle points contribute to the same order. From 
this point of view it seems far from obvious that amplitudes in the mixed regime can be 
obtained by simple analytic continuation and that they contain a single exponential term. 
This issue certainly deserves further investigation. 

There are a few more or less obvious extensions of this work. It seems worthwhile 
exploring the relation between the analytic structure of scattering amplitudes and excited 
states of the auxiliary quantum system for more general cases, including other regimes and 
for a larger number of external gluons. In this context it might be useful to work with an 
alternative version of the non-linear integrable equations [32] in which only the geometric 
cross ratios appear as fundamental parameters. 4 Another question addresses the relation 
between excited states and analytic continuation. For other quantum integrable systems it 
is well-known that there exist excited states that are not related to analytic continuation 
in the parameters. The associated integral equations are of the same form, but without 
any constraints in the numbers N of driving terms. One may wonder whether such excited 
states could also have some interpretation in the context of scattering amplitudes. 

Acknowledgements: We thank Lev Lipatov, Jorg Teschner and Pedro Vieira for inter- 
esting discussions and comments on the manuscript. VS is grateful to the String Theory 
group at the UBC Vancouver and to the Pacific Institute of Theoretical Physics for their 
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A Kinematics in the Regge limit 

In this appendix we shall analyze the behavior of the conformal cross ratios m in the 
ordinary multi-Regge limit 

S » Sl,S2, S3 » — ti,— £2 ~ 1. (A.l) 

When written in light-cone coordinates, the asymptotical behavior of our six momentum 
vectors p% = (pf,p~,pj ,pf) may be parametrized by 

Pi ~ (-ur 1 ,-^,...) 
p 3 ~ (a;, a;" 1 ,...) , 
p 5 ~ (LJ~ b ,LJ b , ...) , 

where w is a parameter that becomes large in the Regge limit and we have only displayed 
the leading terms in the limit oj — > 00. The two exponents a, b are free finite parameters. 

4 Since the contribution A^ liod is most easily expressed in terms of the variable m, however, one must 
still determine this parameter as a function of the cross ratios Ui. 



p 2 ~ (-a;,-^- 1 ,...), 
p 4 ~ (u a ,UJ- a ,. . .) , 
p e ~ ...), 
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Note that all transverse momenta must be finite to satisfy the on-shell condition. The 
dependence of kinematic invariants on oj is now easy to compute 

S ~ LO 2 , S 45 6 ~ W (1+a) , S345 ~ U {l+b) , 

Si ~ 0j( 1 - a 1 , S 2 ~ W (a+6) , S 3 ~ UJ^-V , 

— 1\ ~ W° , — t 2 ~ W° , — S234 ~ (-0° ■ 

Depending on our choice of a and b we can realize various asymptotics of the kinematic 
invariants. If a = b = 1, for example, one finds s 2 ~ s. A regime with s% ~ s 2 ~ S3 is 
reached in case we set a = b = 1/3. Setting a = b = 0, on the other hand, allows us to 
perform a quasi multi-Regge limit in which s 2 ~ 1. 

From the scaling properties of kinematic invariants we can easily deduce the following 
behavior of the cross-ratios (2.4), 

ui ~ w-( a+b ) ~ S2 1 , n 2 ~ w~ (a+fe) ~ S2 1 • (A.2) 

Thereby we have shown that u\ and u 2 vanish at the same order when we take our Regge 
limit. Their ratio u\/u 2 remains constant. It remains to study the properties of 7/3. Using 
the Gram determinant relation to express s in terms of the other eight kinematic invariants 
we find 

^3 = — — — = 1 - f(suk,Si,i+i,i+2) ■ ( A -3) 
S456S345 

After substitution of the gluon momenta we obtain 

/(si,ii, SM+M+2 ) ~ co- {a+b) ~ S2 1 . (A.4) 

Hence, both ui/(l — U3) and u 2 /(l — U3) remain finite for oj — > 00. We have used these two 
ratios to parametrize amplitudes in the Regge limit. 

B Hirota equations 

There exists another important set of functions T s ^ m that appear in the context of gluon 
scattering. Our discussion in this appendix will be kept general, i.e. it applies to an arbitrary 
number n of gluons. The functions T s>m = T S)m {6) are parametrized by two integers s,m. 
They are subject to the following set of Hirota equations, 

-^s,m-^4— s m — -^4— s,m+l-^s,m— 1 -^s+l,m-^s— l,m > Q^"^-) 
for s = 1,2,3. The superscipts ± indicate shifts of the spectral parameter, i.e. 

/±(0) = f\±\e) with /H( ) = f(9 + iml) . 
We complete the description of the functions T SyTn by a list of boundary conditions. 

T s _i = T Sjn _4 = 0, where s = l,2,3 (B.2) 
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and 

T Si o = To >m = T^ m = 1 . (B.3) 

Here, n is the number of gluons. After imposing boundary conditions, we are left with 
3(n — 5) functions. 

The T s ^ m are related to y s m we have used throughout our discussion through the 
relations 

v J- s,m+l-L 4-s,m-l /tj ,i\ 

s+l,m^ s— l.m 

Inserting these into the Hirota equations (B.l) we obtain 

^SjTnXi— s,m (1 + Xs,m+l)(l + Y— s,m— l) 



(1 + is+l,m)(l + Ys-1, 



(B.5) 



where s = 1, 2, 3 and m = 1, 2, . . . , n — 5. The behavior of T S)m on the boundary translate 
into the following boundary conditions for Y s>m , 

Yo t m = y 4jm = oo, and Y sfi = Y Sjn _ 4 = with s = 1,2,3. (B.6) 

Using the set of Y— functions one may relate the parameters m m , cp m , C m of the Y— system 
through cross ratios. In order to do so, one has to invert the following relations 

„2 2 / 1 \ / T+ T 

p-k-l,p+k- L p-k-2,p+k~l _ ( _|_ 1 I _ / ■ L 2,2k-1 2,2k-l 



2 2 \ ~ ' V I IT 1 ! -1 

X p-k-2,p+k X p-k-l,p+k-l V *2,2k-l / 0=i^(2p-l) \ J 2,2fcJ2,2fc-3 y 6)=i2L ( 2p _i) 

(B.7) 

Our equations (3.8) are recovered as a special case when n = 6 with k = 1 and p = —1, 0, 1. 
We have pointed out before that the Y-functions satisfy the usual non-linear equations 
as long as their argument 9 remains inside the strip |lm(0)| < ir/4. Whenever 9 crosses 
one of the lines Im(9) = ir7r/4, one has to add contributions from poles of the kernel 
functions K^. To avoid these complications, we can also solve the original equations inside 
the fundamental strip and extend the solutions to the complex plane recursively, with the 
help of the Y — system equation, i.e. 

y (1 + ^s,m+l)(l + ^4-s,m-l) /g gN 

~ Yll im (l + 1/Y} l \j(l + 1/Y}% m ) ' 

We have used this recursive construction of Y s = Y s \ is section 3.1 to obtain the formula 
(3.9) for the cross ratio M3. 

C Comments on numerical computations 

In this work, numerical calculations were used to find the path of the analytical continuation 
in the space of parameters m,(j),C and to check which solutions of Y = —1 cross the 
integration contour. Most of this was done with the help of Mat hematic a. 

In order to determine the Y — functions for fixed m, (j) and C we use an algorithm similar 
to one presented in [29] for the case of AdS^. It involves solving the non-linear integral 
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equations (3.5) and (3.6) perturbatively for large enough m. The leading solution is given 
by expressions that do not contain integrals, namely our eqs. (3.12) and (3.13). To obtain 
the solution with better precision we insert the k th iteration into the integrals to obtain, 
i.e. schematically 

log>f +1) (#) = -m 4 cosh(0-i ( / ) ) + C 4 + (K ii *log(l + y i (fc) ))(e). (C.l) 

We perform k such iterations until the Y— functions stabilize. As a result we obtain 
Y— functions of real 6, or rather a discrete subset thereof. Moreover, using the integral 
equations we can calculate the Y— functions for complex value of 6 by convoluting the 
kernel functions -fQ with the Y— functions of real 8. Due to existence of the poles in the 
kernels, we must limit this evaluation of the ^-functions in the complex 0-plane to the strip 
| Im#| < For values of 9 outside this strip, we recursively use Hirota equations (B.8). In 
this way, we can construct Y — functions on the whole complex 0— plane. 

Next, we evaluate the Y— functions for special values of 6 = ±7ri/4, 37ri/4 which de- 
termine the cross ratios Uj defined in eq. (3.8). Actually, for these values of 9 the numer- 
ical algorithms are not well convergent and it is easier at these points to interpolate the 
Y— functions. But at the end one can relate values of m, (ft and C to values U{. In order to 
find Y — functions along our analytical continuations contour defined by eq. (2.11) one has to 
invert the relation between the parameters m, (ft an d C and the cross ratios U{. The values 
of C along the path can be found analytically, namely it is given by eq. (5.3). The other 
two parameters m and (ft, on the other hand, require numerical calculations. An example 
of such path is shown in figs. 2,3. Here, it has been performed by solving a set of nonlinear 
equations (3.8) using the Newton method. 

Furthermore, one can also determine the positions of solutions to the equations Y^ = — 1 
numerically, see figs. 4,5. These solutions are moving while we are continuing along the 
contour in parameter space. In the process of this continuation, some solutions of I3 = — 1 
and Y\ = — 1 cross the real line. This happens roughly when ip ~ ir/2. At this point we 
must replace the original non-linear equations by those spelled out in eqs. (5.4) and (5.5). 

The numerical results show that during analytical continuations only from of the so- 
lutions to Yi = — 1 cross the integral contour, see figs. 4 and 5. Moreover, to a very good 
approximation, the imaginary parts of m and (ft behave linearly with tp while their real parts 
change quadratically, compare also fig. 2. At the end of the contour path, m and (ft go to 
values which are captured by our eqs. (5.13) and (5.14). 
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